Draft version April 18, 2013 

Preprint typeset using L^-T^X style cmulatcapj v. 12/16/11 



REDUCING SYSTEMATIC ERROR IN CLUSTER SCALE WEAK LENSING 

Yousuke Utsumi 1,2 , Satoshi Miyazaki 1,3 , Margaret J. Geller 4 , Ian P. Dell' Antonio 5 , Masamune Oguri 6 , Michael 

J. Kurtz 4 , Takashi Hamana 1 , and Daniel G. Fabricant 4 
Draft version April 18, 2013 

ABSTRACT 

Weak lensing provides an important route toward collecting samples of clusters of galaxies selected 
by mass. Subtle systematic errors in image reduction can compromise the power of this technique. We 
use the B-mode signal to quantify this systematic error and to test methods for reducing this error. We 
show that two procedures are efficient in suppressing systematic error in the B-mode: (1) refinement 
of the mosaic CCD warping procedure to conform to absolute celestial coordinates and (2) truncation 
of the smoothing procedure on a scale of 10'. Application of these procedures reduces the systematic 
error to 20% of its original amplitude. We provide an analytic expression for the distribution of the 
highest peaks in noise maps that can be used to estimate the fraction of false peaks in the weak lensing 
/c-S/N maps as a function of the detection threshold. Based on this analysis we select a threshold S/N 
= 4.56 for identifying an uncontaminated set of weak lensing peaks in two test fields covering a total 
area of ~ 3deg 2 . Taken together these fields contain seven peaks above the threshold. Among these, 
six are probable systems of galaxies and one is a superposition. We confirm the reliability of these 
peaks with dense redshift surveys, x-ray and imaging observations. The systematic error reduction 
procedures we apply are general and can be applied to future large-area weak lensing surveys. Our 
high peak analysis suggests that with a S/N threshold of 4.5, there should be only 2.7 spurious weak 
lensing peaks even in an area of 1000 deg 2 where we expect ~ 2000 peaks based on our Subaru fields. 



o 

(N 

< 



43 

6 

5h , 1. INTRODUCTION 
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Weak lensing is a fundamental tool of modern cosmology. A weak-lensing map provides a weighted "picture" of the 

projected surface mass density and thus a rou t e to identifying clusters of galaxies selected by mass Ce.g.lMivazaki et all 

2002a; D ahle et al.ll2003t iSchirmer et al.ll200l iHetterscheidt et aLll2005t IWittman et al.ll2006t iGavazzi fc Soucailll2007t 

Miyazaki et al. 2007; S han et aLll2012| ). Subtle issues limit the application of weak lensing maps as sources of cluster 
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catalogs. A serious astrophysical limitation is the projection of large-scale structure along the line-of-sight (see e.g 
^D ' iWhite et~al . 2002; Ham ana et al. 2004, 2012). An observational limitation is the presence of systematic errors in the 



I/"") ' maps dGeller et al 



2010; Kurt z et al. 20121) . Here w e focus on a method of reducing these errors. 



VO , IWittman et al.l ( 2001f) and iMivazaki et all (|2002aT ) first identified rich clusters from weak lensing maps. Since then 



I larger and larger weak lensing surveys with increasingly sophisticated reduct ion techniques have led to a steep increase 



in the number of clusters initially detected in weak lensin 



Schirmer et~aLll2007t IGavazzi fc Soucailll2007l IMivazaki et al 



g maps (Hcttcrschcidt et al. 2005; Wittman et al. 200 



20071; iBcllaga mba etal.ll2011l; IShan et al.ll2012r i. iShan et al.l ([20121 ) used the CFHT to carry out the largest area survey 



120071: iMaturi et al.ll2 007: Be rge et al.ll2008l : lDiet rich ct al. 



to date covering 64 deg . Among the 301 weak lensing peaks with signal-to- noise ratio greater than 3.5, only 126 have 
a corresponding brightest cluster galaxy within 3' of the weak lensing peak. Shan ct al. (2012) conclude that many of 
the weak lensing peaks, even above this threshold, ar e probably noise . 

In a detailed reanalysis of the 2.1deg 2 GTO field of IMivazaki et al.l (j2002all2"o"07l ); Kurtz et all (|2012l ) show that the 
highest peak of the B-mode map (in this map the sources ellipticities are rotated by 45° and there should be no signal) 
C3 • has signal-to-nois e, S/N ~ 4.25 and there are four peaks in the B-mode map with S/N> 3.7. It is interesting that 
IShan et al.1 (120121 ) find a si milar numbe r of high-significance peaks per unit area in their B-mode maps. Analysis of 
the B-mode map prompted Kurtz et al. (2012) to set their threshold for weak lensing cluster detection at 4.25. Above 
this threshold 2/3 of the weak lensing peaks correspond to individual massive clusters. 

Here we investigate possible underlying sources of the high peaks in the B-mode maps. Our goal is reduction of 
the amplitude and frequency of these peaks as a basis for construction of more robust catalogs of weak lensing cluster 
detections. We apply our r evised procedures to two fields observed with the Subaru Telescope: the GTO field of 
Miyazaki ct al. (2002a], I2007D and a small portion of the DLS F2 field reobserved with Subaru. 

In Section 2, we review the standard analysis of the Subaru imaging data. In Section 3 we review the construction of 
the galaxy catalog. In Section 4 we describe the basic weak lensing analysis we use. In Section 5 we develop a revised 
procedure for reducing systematics in the B-mode map. We demonstrate that there is excess power in the B-mode 
map on small and large scales relative to noise maps. Two procedures are effective in suppressing systematic in the 
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TABLE 1 
Observation Log. 



Name 


Ra 


Dec 


Exp. time (sec) 


Observation date 


Filter 


Seeing 


Limit mag 


GTO.0 


16:04:44.5 


+43:11:12 


2100 


2001/4/23, 25 


Rc 


0.66 


25.38 


GTO.l 


16:07:39.0 


+43:12:19 


2400 


2001/4/23, 25 


Rc 


0.77 


25.28 


GTO.2 


16:07:39.0 


+43:37:37 


1800 


2001/4/23 


Rc 


0.72 


25.44 


GTO.3 


16:04:44.6 


+43:36:30 


1800 


2001/4/24 


Rc 


0.71 


25.24 


GTO.4 


16:01:46.8 


+43:37:37 


2250 


2001/4/24, 25, 5/18 


Rc 


0.65 


25.17 


GTO.5 


16:01:46.8 


+43:12:19 


2250 


2001/4/24, 25 


Rc 


0.70 


24.90 


GTO.6 


16:01:46.8 


+42:47:01 


1800 


2001/4/24, 25 


Rc 


0.66 


25.19 


GTO.7 


16:04:43.0 


+42:47:01 


1800 


2001/4/25 


Rc 


0.65 


25.18 


GTO.8 


16:07:39.0 


+42:47:01 


1800 


2001/4/25 


Rc 


0.64 


25.32 


DLSF2J0 


09:16:30.9 


+29:17:40 


3600 


2008/1/8 


i' 


0.79 


25.50 


DLSF2J1 


09:16:31.0 


+29:39:40 


3600 


2008/1/8 


i' 


0.72 


25.52 


DLSF2J2 


09:16:31.0 


+30:01:40 


3600 


2008/1/8 


%' 


0.72 


25.33 


DLSF2J3 


09:16:31.0 


+30:23:41 


3600 


2008/1/8 


i> 


0.76 


25.51 



a The GTO pointings fill a 3 X 3 square. They are numbered so that the central field is and numbers go from 1 in the central left (east) 
field clockwise to 8 in the southeast corner. The DLS pointings fill a 1 X 4 rectangle: is the southernmost pointing and the numbers 
increase to the north. For the GTO field the numbering is the same as in Kurtz ct al. (2012) 

b Detail of the limiting magnitude as described in Section [3] 

B-mode map: (1) warping the image onto absolute sky coordinates and (2) introducing a large-scale (10') cutoff in 
the smoothing kernel applied to construct the weak lensing map. We define the weak lensing peaks and investigate 
thresholds in Section 6. In Section 7 we note some differences in our analysis of the GTO and DLS Subaru observations. 
Section 8 provides an analytic expression for the probability of finding the highest peak in the noise maps as spurious 
peaks in a weak lensing map. In Section 9 we test the lensing maps of our two fields by comparing the peaks with 
systems identified in dense redshift surveys. We demonstrate that a threshold S/N = 4.56 produces a robust catalog 
of peaks with few false positives. 

Unless otherwise stated, we adopt the standard concordance cosmology (h = 0.72, Q m = 0.27, tt\ = 0.73) and the 
AB magnitude system throughout this paper. 

2. IMAGING DATA & IMAGE REDUCTION 

In this section we describe a revised reduction procedure that reduces systematic errors in weak lensing ma ps. As ex- 
ample s of the revised procedure, we use two regions imaged with Suprime-Cam on the Suba ru Telescope ( Miyazak i et al.l 
l2002bD . The GTO field covers ~ 2 deg 2 , centered on (16:04, +43:12) and was observed bv lMivazaki et a l. (2002 a]). We 
re-ob served the western ~ 1 deg 2 portion, (9:16, +30:00), of the Deep Lens Survey F2 field (DLS F2 [Wittma n et all 
l2002t) . The new observations were taken in 0.75" seeing. The GTO field consists of 9 Suprime-Cam pointings; the 
DLS observations consist of 4 Suprime-Cam pointings. Table Q] shows the numbering scheme for Subaru pointings in 
both the GTO and DLS fields. 

We retrieved data from the SMOKA data archive server (jBaba et al.l I2002D . These data satisfy the condition 
FILTER=W-C-RC and PSF_SIGMA < 0.9 to obtain good i mage quality. PSF SIGM A is the typical FWHM (arcsec) 
of stellar images evaluated by SMOKA using SExtractor ()Bertin fc Arnouti [1996) • If the total exposure exceeds 
2,000 seconds, we select frames with the smallest PSF_SIGMA and reject frames with larger PSF_SIGMA even if 
PSF.SIGMA is < 0.9". As a result, we reject about 6 of 11, 10 of 15 and 1 of 5 exposures for GTO.0, GTO.4, and 
GTO_5, respectively. Table [1] gives the total exposure times, observation dates, filters, typical seeing, and effective 
limiting magnitude for each subfield. Imaging observation for the DLS F2 field were made with Suprime-Cam in 
January 1, 2008. Each of the 4 subfields is a 3600 sec exposures taken under stable seeing conditions. Table [1] lists the 
FWHM of the PSF on stacked image s in each of the f our subfields . 

The basic image reduction follows Miyazak i et al.l (|2002al I2007D and iKurtz et al.l (|2012|) . However, we make one 
major revision: we use independent astrometric data to register the images accurately with respect to the celestial 
coordinate while avoiding the introduction of any additional image warp. 

The preliminary reduction of the raw image on each CCDs is standard. We use SExtractor to identify objects on the 
overscan-subtracted, flat-fielded, sky-subtracted image. We identify stellar objects from their easily identifiable linear 
size-magnitude relation. We manually select the stellar object catalog using this relation. 

After completing the preliminary reduction, there are three additional steps in obtaining the final mosaic of stacked 
images; one of the steps is new. First we determine the accurate locati on of each CCD in instrume nt coordinates 
(Section l2~Tj) . Second, we apply the image alignment warping procedure of iMiyazaki et al.l (|2002al , 120071 ) ( Section |2~2|) . 
Finally we introduce a new additional warp "Registration onto the celestial coordinate" ( Section 12. 3|) . 

2.1. Solving the Basic Mosaic Geometry 

We parameterize the accurate location of each CCD in instrument coordinates by the displacement and the rotation 
of each CCD (Ax, Ay, Acj)) c , the telescope pointing offset between dithered exposures (AA^, Ay, A$) e , and the optical 
distortion of the wide field corrector lens. The distortion parameters are well modeled by a 4th-order polynomial 



function: 



R — r 

= ar + br 2 + cr 3 + dr A (1) 



where R and r are distances from the optical axis in units of pixels on the CCD. R is the original instrument coordinate 
and r is the distortion-free coordinate. 

To construct a mosaic we need a flux calibration to correct for the differences in sensitivity from chip to chip (A/ c ). 
These differences originate from differences in amplifier gain, differences in the normalization of flat fielding, and 
from corrections for the transmission of the atmosphere from exposure to exposure A/ e . We derive this calibration 
by comparing the fluxes for unsaturated stars in overlapping regions of the frames. We call these parameters the 
"mosaicking rule" . 

Operationally, we reduce the rotational component (A<p, A$) to the linearized form of the rotational matrix consisting 
of four components ( (4>ij), (*&ij)) where i,j £ (1,2). Because this matrix can treat not only rotation but also expansion, 
we include the first coefficient of the distortion parameter in this matrix. Varying the basic mosaic geometry parameters 
(Ace, Ay, 4>ij,Af) c , (AX, AY, <%, A/) e , (b, c, d), we minimize 



a 2 = v y ^ -ui i + wi -h > (2) 




where iff and f\ are the e-th exposure's position in distortion-free coordinates and the corrected flux for i-th unsatu- 
rated star and o~i, i € (rj, /) are the typical error in the measurements of position and flux, respectively. We fix param- 
eters for the first exposure of the series and the left bottom chip, i.e. (Ax, Ay, A(j>, A/)o = 0, (AX, AY, A$, A/)o = 
as constraints. 

This mosaicking rule is calculated by imcat dKaiser et al.l ll995). The best fit parameters are obtained by minimizing 
the positional difference of unsaturated control stars (70 ~ 100 stars per CCD) for each exposure. 

The residual alignment error in this procedure is ~ 0.5 pixel rms (0.1"), still large compared with typical distortions 
of the sources in weak lensing. Many stacking pipelines adopt this "typical stacking" warp only. We note that a typical 
stellar size is 0.6"; it should thus be possible to align the stellar positions to < 0.06", one tenth of the stellar size. The 
accuracy of the "typical stacking" is insufficient and probably below the level that can be obtained from the data. 

2.2. Jelly Focal Plane Warp 

Residual positional differences from exposure to exposure in the typical stacking procedure in Section [2. II may result 
from non-symmetric features of optical distortion (possibly due to imperfect optical alignment) not considered in the 
model. Differential atmospheric dispersion is another source of asymmetry. 

We parameterize the residuals relative to one another by a polynomial function of the field position; we expect the 
residuals to be continuous over the field of view. The residuals are 

6 I 6 1 

AC = Yl E ai m , e e- m V m , A v = yy h m ^ l - m rr (3) 

1=0 m=0 1=0 m=0 

where ff = (£,T)). We then obtain the coefficients ai m , e and 6; m ,e by minimizing the variance of the residuals. Each 
individual CCD image is "warped" usi ng this polynom i al correction prio r to the stacking. This process reduces the 
alignmen t error to ~ . 05 pix el (0.014"). IMivazaki et al.1 (j2002aLl2007h and lKurtz et all (|2012f l used this transformation 
following Kai ser et al.l (|1998t ). We refer to this trans formation as "the original wa rp" , and we call the res ults of this 
warp "the original stacking". The result is similar to IMivazaki et al.l ()2002al [2007T l and lKurtz et al.l (|2012D . 

2.3. Registration onto the Celestial Coordinate 

The third transformation is new and we apply it as the final step in the reduction. The result of Section 12.21 is 
that all exposures except for the first one are registered on the first exposure with an accuracy of ~ 0.05 pixel, the 
tolerance we seek. However, there is no guarantee that the geometry of the first exposure is well matched to the 
celestial coordinate. 

In practice, the absolute rotation of the instrument rotator is less well-known than its rotational angle relative to 
its zero point. There is no mechanical system to determine the zero point rotational angle for the rotator accurately 
for Suprime-Cam; the rotational angle relative to the origin is, however, very well known. The possibility of this zero 
point offset suggests that the final Suprime-Cam image will be rotated relative to the absolute celestial coordinate. 

To align the final image with the celestial co ordinate fram e (not some artificial distorted frame, but a physical 
frame), we use an external catalog, USNO-A2.0 (|Monetlll998t ). This astrometric standard catalog has an accuracy of 
0.4 arcsec although USNO-A2.0 is not deep enough to match the Suprime-Cam catalogs. The overlapping magnitude 
range includes moderately bright stars; these stars are saturated but do not have blooming pixels. It is thus difficult to 
measure their total intensity, but it is possible to measure their position because they have circular symmetry. We do 
not use the SDSS star catalog because it does not cover the enough of the sky to allow us to to extend our procedure 
outside the GTO field. We confirm that using the SDSS external astrometric catalog makes no significant difference 
in the warp relative to use of the USNO-A2.0. 



In the same manner as the transformation in Section 12.21 we parameterise the residuals (Ace, Ay) of bright stars 
(they are slightly separated from the linear stellar sequence but have a comparable positional determination to the 
unsaturated stars provided there are no blooming pixels) in the first exposure with respect to the USNO-A2.0 celestial 
coordinate as a polynomial function of field position x — (x. y) (Equation|3|). This procedure yields absolute astrometric 
object positions accurate to ~ 0.2". The internal relative alignment of stars remains at the level of ~ 0.05 pixel. We 
refer to this transformation as "the new warp" and we call the result of this warp "the new stacking" . 

We present a schematic picture to describe these three steps of the image reduction in Figure Q] and summarize the 
three step image reduction pipeline (Section 12. 1[ I2.2[) compared to original one and typical one in Table [2] Once 
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Fig. 1. — Schematic of the three steps in the mosaic warp 



TABLE 2 
Summary of the warp rules 





New 


Original 


Typical 


Registration 


Distortion f(R) 

Jelly-warp 

Align to sky 


yes 
yes 
yes 


yes 
yes 
no 


yes 
no 
no 


~ 0.5 pixel 
~ 0.05 pixel 
~ 0.05 pixel 


Astrometry 


~ 0.2 arcsec 


uncontrolled 


depending on pipelines 





we have all of the parameters that describe these three steps, we carry out the image warp all at once rather than 
applying the transformation step by step. This integration of the procedures is necessary to minimize the artifacts 
due to the interpolation of pixels. We estimate pixel values via 3rd order bi-polynomial interpolation of 4 x 4 source 
pixels. We employ the arithmetic mean for stacking. 



3. GALAXY CATALOG 

To construct the galaxy catalog from the mosaic image, we use the imcat tools "hfindpeaks" (object detection) and 
u apphot" (photometry). We use the half-light radius "?V as the size of an object. Using this size, we can distinguish 
galaxies from stars by requiring that rn for galaxies satisfy the relation: 

r h >r* h + a rl , (4) 

where r£ and <r r * are the half-light radius of a a stellar image and its rms, respectively. 

We merge our catalog of stars with the SDSS star catalog to calibrate the photometric zeropoint (see. Equations 
[5]). Typically, the number density of stars to the depth of SDSS is several arcmin -2 ; the precision of the astrometric 
catalog is < 1". 

We determ i ne th e photometric zeropoints for each subfield from the the SDSS photometric catalog following 
lUtsumi et al.l ([2010D . The transformation equation between the photometric systems is: 

Rcsubam - r' SDSS = -0.314(r^ DSS - i' SD ss) ~ °- 011 for r SDSS ~ i'sDSS < °- 5 

*S«6aru ~ *S£>ss — 0. 102(i s . £)S . s . — r SDSS ) -(- 0.002 for - 0.5 < i SDSS ~ r SDSS (5) 

To assess the depth of our image s, we follow the standard procedure of SDFRED to derive the limiting magnitude 
()Yagi et alj|2002t lOuchi et~aLl l2004'). We give the 5er limiting magnitude in a 2" photometric aperture in Table [T] 

To remove low signal-to-noise regions around sufficiently bright stars, we need a procedure for constructing masks. 
Stars brighter than 20 th magnitude are usually saturated in a Suprime-Cam typical exposure with ~ 100 sec inte- 
gration. We measure the actual PSF of Suprime-Cam over a wide dynamic range by compiling the PSF of faint and 
bright stellar objects. We divide the stellar samples into subsamples segregated by R-band magnitude, faint (20-20.5 
magnitude), intermediate (18-18.5 magnitude) and bright (13.75-14.25 magnitude), combine them by taking the arith- 
metic mean, and then measure their radial profile. Although we take a mean profile of bright stars, we also reject 
exceedingly bright pixels by choosing an appropriate threshold for each images in order to avoid saturated areas or 
blooming regions. 

For stars with R < 14, there are not enough stars to determine the mean PSF. Thus we consider an optical model 
of Suprime-Cam to determine the size of the surrounding halo. The halo surrounding a bright star is a reflection from 
optical glass components in front of the CCDs. The nearest and second nearest optical elements close to the CCDs 
consist of the entrance window to the dewar and the filter; their distances from the CCDs are 10.0mm and 34.5mm, 
respectively. For both elements, the thickness and effective reflection index is 15.0mm and 1.460280, respectively. The 
focal ratio Suprime-Cam is 1.86. Thus light reflected on the first surface, the backside of the dewar window, should 
be located 71.3 arcsec away from original position. Light from the other surfaces appears 144 , 212, 284 arcsec away 
from the original position. Q This model is consiste nt w ith the data and we use it to derive the mask radii. Adopted 
mask radii for brighter stars are shown in Appendix IA11 

We mask some stars with R < 11 by visual inspection to remove obvious artifacts. Appendix IA2I and IA3I show the 
masks. 

4. BASIC WEAK LENSING ANALYSIS 

To construct an initial weak lensing map we follow the procedures described in iMivazaki et al.l ()200 2a. 2007]) and 
iKurtz et al.l ((2012). In this section we review the selection of sources and the construction of the mass, noise, and 
B-mode maps. In discussing the B-mode map, we motivate the investigation (Section [5]) of sources of systematic error 
in this map. 

4.1. Sources 

We use 23 < m < 26 galaxies with detection significance v > 10 as source galaxies. We eliminate high ellipticity 
objects with |e| > 0.6 because these objects are sometimes blended. This magnitude range gives an average magnitude 
of 24.5 for the sources in both the GTO and DLS fields. The average source density is 34.5 arcmin -2 and 31.6 arcmin -2 
for the GTO and DLS fields, respectively. Figure [2] shows number density maps of sources galaxies for weak lensing 
analysis. 

Bec ause we do not have measured redshifts for the source galaxies, we use the fitting formula from lSchrabback et al.l 
(J2010D for galaxy redshifts as a function of magnitude in the range 23 < isi4 < 25. With this relation we estimate an 
average redshift of 1.2 for the DLS field where we have i-band data. For the GTO field where we have only .R-band 
date, there is no suitable redshift-magnitude relation. We thus assume the same average redshift as the DLS field. If 
the actual redshift is 1.0 rather than assumed 1.2, we introduce an ~ 10% error in the estimate of the distance ratio. 
We note that this error in mass estimation is only important for absolute calibration. 

To derive the ellipticities of the sources, we use getshapes, efit, and ecorrect in imcat. We evaluate the re- 
circularization correction using stars in the field of view and we model the PSF dis tortion as a fifth order bi-polynomial 
as a function of position. The P 1 correction follows KSB+ (|Hoekstra et al.|[T998D who adopt a galaxy size dependent 
correction. The actual 7-field is discretized and noisy because it can only be estimated at the positions of galaxies. 
Note that we do not weight the estimated shear value according to galaxy's magnitude. We apply a gaussian smoothing 

7 Mivazaki ct al. (2002b) says the distance between the window and the filter is 14.5 mm but this turns out to be modified to be 9.5 mm 
in the actual implementation of Suprime-Cam. The translation of the flat glass (filter), however, does not affect any optical performance. 
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Fig. 2. — Surface number density of weak lensing sources. For The GTO field (left) the sources have 23 < R c < 26; for the DLS field 
(right) the sources are in the range 23 < i' < 26. 



filter on the 7 fiel d following lOguri fc Hamanal (|2011h . 
IKurtzet al l (pOll . 



At first, we adopt a smoothing length, 9q — 1.5', following 



4.2. Constructing the k-S/N Map 

We can now construct the k-S/N map to search for mass concentrations based on the significance of the lensing 
signal. We derive the k map directly from the 7 map of the sources in Section |4. II We divide this map by the noise 
map to obtain the k-S/N map. 

We use Monte Carlo realizations to construct the noise maps. We rotate the orientation of each source galaxy by a 
random angle. The randomized catalog is then the basis for the noise map. We repeat the Monte Carlo process 100 
times, and then calculate the rms at each grid point of the noise map to estimate the noise at that point. The k-S/N 
map then reflects the noise from intrinsic ellipticity noise and from variation in the number of galaxies. 

In the following sections (Section 14.31 14.4]) . we examine the systematic error in our weak lensing map by looking at 
the S/N distribution of the noise and B-mode maps. We demonstrate that the noise map obeys the expected Gaussian 
S/N distribution; however, the original B-mode map does not. 

4.3. Testing the Noise Map 

In the absence of signal, the S/N distribution for the noise map should match a standard normal distribution which 
has no free parameters. We test this expectation by constructing a "Noise S/N map" for each of the 100 Monte Carlo 
realizations. Pixel values on each of these map are strongly correlated spatially because the maps are smoothed with 
the Gaussian kernel. However a series of pixel values at a fixed position in the 100 "Noise S/N maps" should obey a 
normal distribution. 

We show the S/N histogram for randomized realizations in Figure |3l The random realizations have some scatter 
that we use to evaluate the presence of systematic errors in the B-mode map. To quantify the scatter, we define the 
"Deviation" : 



Deviation 



N 



obs 



N? 



vw 



(6) 



where iVj = dN/d(S/N)i is a number of pixels per d(S/N)i bin in the observed (obs) and in the standard normal 
distribution (norm). The error is the standard deviation. We obtain the merit function as a kind of reduced chi- 

X 



squared, \ 2 > from a sum of the squared deviations dividing by the the number degree of freedom (the number of 
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Fig. 3. — [Top] S/N distributions for 100 noise maps in the GTO field, which are obtained by randomizing orientations of source galaxies. 
Different points indicate different realizations. The dashed line shows the expected normal distribution. [Bottom] panel shows the deviation 
from the normal distribution (Equation [SJ . The sum of these deviations divided by the degree-of- freedom is a measure of the systematic 
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Figure [4] shows the normalized cumulative distribution of \ 2 for all random realizations. This wide range, < \ 2 < 65, 
gives the maximum allowable difference between the noise and B-mode maps. 

4.4. Testing the B-mode Map 

In this section, we quantify the quality of the B-mode map. The standard B-mode map results from rotating each 
shear |7|exp(2z0) by 7r/4, i.e. 71 + 172 — > —72 + iji, where 9 is the position angle for each source. Because weak 
lensing does not produce a B-mode, we can evaluate the level of systematic error inh erent in the mass reconstruction 
procedure by comparing the B-mode map with the noise map ()Schneider et al.ll2002f) . 

Figure [S] shows the pixel distribution for the B-mode maps with the original, typical, and new stacking (Table |2J. 
The S/N distribution has a large deviation relative to the standard normal distribution. The x 2 — 257 (1204) for the 
original (typical) stacking. 

A possible origin of these B-mode artifacts is misalignment during the stacking procedure. Assuming we are stacking 
circular gaussian objects, a 10% positional inaccuracy generates ~ 1% artificial ellipticities. The positional difference 
applies to individual galaxies and their neighbors. The smoothing procedure then leads to detection of these ~ 1% 
ellipticities as a lensing signal. Thus stars must be registered internally with an accuracy of 1% level (^0.05 pixel). 

Next we test procedures for reducing this kind of artificial signal in the B-mode map: (1) SkyFrame, (2) Wide 
FFTBox (3) Masking the edge of the field, (4) Higher order PSF anisotropy correction and (5) Smoothing cut off. 




Fig. 4. — Normalized cumulative distribution of the merit function \ 2 (Equation from the 100 realizations of the noise maps. 



5. REDUCING THE B-MODE SIGNAL 

Systemic issues which might produce an artificial signal in the B-mode map include errors in the physical coordinate 
system, excess large-scale power, boundary effects, and errors in the PSF anisotropy correction. 

Sky Frame — If the stacked image is bent with respect to the physical coordinate system, the E-mode signal contami- 
nates the B-mode. The tilt of the coordinates introduces a rotation of the sources. 

The "new stacking" ( Section 12. 3[) eliminates this bending. The "new stacking" image is aligned relative to the sky 
coordinate within ~ 1 pixel; the internal registration accuracy is ~ 0.05 pixel. 

Wide FFTbox — There are two ways to perform the KSB analysis; (1) convolve the shear field with the complex kernel 
or (2) perform the convolution in the Fourier domain. Mathematically these procedures are identical. However, in 
practice for operation i n the Fourier dom ain, boundary effects can be a problem. To avoid boundary effects present in 
the original reduction (jKurtz et al.H2012T ). we embedded the shear signal for the E/B-mode maps within a calculation 
box twice the size of the field with zero padding outside the map area. 

Masking t he edge of the field — A t distances more than 15' from the optical axis the field of view of Suprime-Cam is 
vignetted (Miyazaki ct al. 2002b). The distortion is also large, > 10", at the edge of the field. We thus try masking 
the region outside 15'. 

Smoothing cut off — To examine the origin of the remaining B-mode signal, we examine the power spectrum of the 
B-mode. 

There is a significant excess in the power spectrum of the B-mode map relative to the noise map on large scales 
(> 10'). This excess is comparable with the power in the E-mode (Figure EJ. The origin of this signal in the B-mode 
is unclear, but it should not be present. 




Fig. 5. — [Top] Histogram of S/N values for B-mode maps derived from the new result (circles) as well as original (triangles) and typical 
(squares) stacking. The dashed line shows the expected normal distribution if there is no systematic error. [Bottom] Same as [3] but for the 
three kinds of B-mode. 
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To suppress this large scale systematic error in the B-mode map, we truncate the smoothing kernel at a scale of 10'. 
Figure [5] shows the result. The truncation reduces the B-mode signal. 



Fourier modes of S/N map 



le+09 F 



le+08 






le+07 : 



le+06 r 



- 1 1 r - 



- 1 1 1 1 — i — i — r - 



3e+08 



Typical separation of stars 



5e+07 




ft, l«fe^A»4-"' 



Minor width of CCD 

Original 
Sky Frame / WideFFTBox 

Noise 
Noise 



0.1 



1 



10 



Wavelength of the fourier mode [arcmin] 



Fig. 6. — Power spectrum of the noise maps (short dashed and dotted lines) and B-mode maps constructed with various procedures (all 
other curves as noted in the figure). The power spectrum for the B-mode has an excess relative to the noise maps on smaller scale less 
than the typical separation of stars and on scales larger than 10 '. For scales smaller than 1.5 ', smoothing decreases the amplitude of the 
power spectrum. The inset shows the power spectrum in the range 1.5 to 3 '. The SkyFrame procedure minimizes the excess in this range. 
A higher order PSF correction has no impact. For scales of 1 - 10 ', there is no significant excess of in the B-mode power spectrum. For 
scales larger than 10', all procedures fail to reduce the B-mode excess. Because the typical scale of a cluster of galaxies is few arcminutcs, 
these scale can safely be filtered out. 



Higher Order PSF anisotropy correction — To correct the PSF anisotropy induced by atmospheric disturbance or optical 
aberration, we applied a fifth order polynomial PSF model as a function of position. If the order of the polynomial is 
not sufficient to remove the systematic variations in the PSF, a B-mode is produced. To investigate this issue, we test 
a sixth-order polynomial. 

5.1. Summary of the Tests 

Table [3] shows the result of these procedures. The "SkyFrame" registration reduces the B-mode power on small 
scales and the "Smoothing Cut-off" reduces the power on large scales (Figure |BJ). These two procedures have the 
largest effect on reduction of systematics in the B-mode map. 

The "WideFFTBox" also reduces the B-mode but not significantly. We thus use direct convolution instead of the 
FFT implementation. 

"Masking the edge of the field" and "Higher order PSF anisotropy correction" have little effect for the GTO field. 
"Masking the edge of the field" does not suppress systematics in the B-mode map. The same is true of the higher 
order corrections for PSF anisotropy. 

Based on these results, we include 

• SkyFrame image 

• Smoothing Cut-off at 10 arcmin 
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TABLE 3 
Summary of the Tests 



Tests 


x 2 


Typical stacking image 






1204.7 


Original result^ 






257.7 


SkyFame 






90.7 


SkyFrame+Wide FFTBox 






81.9 


SkyFrame+Wide FFTBox+6th 






84.3 


SkyFrame+Wide FFTBox+Masking 


the ed 


gc 


86.3 


SkyFramc+Dircct Convolution+Smoothing 


cut off 


54.4 



IKurtz et~aTI 112013) 

b Result values from the merit function defined as Equation [7] according to adopted tests. 

• Direct convolution instead of using Wide FFT Box 

in the construction of the maps. We show the result of this procedure in Figure [5] The pixel distribution of the B-mode 
S/N map indicates that residual systematics are reduced to a reasonably low level if we adopt appropriate procedures. 
To assess whether the deviation of the B-mode signal is consistent with random realizations, we compare x 2 values 
with the distribution for each random realization. The maximum x 2 is larger than what we derive from the B-mode 
signal for the GTO field (Figure 0}, x 2 = 54.4, much less than the original 257. Thus we confirm that the our revised 
lensing map is substantially less affected by systematic error. Figure [7] shows the resulting mass map with the "new 
warp" and Figure [8] shows the resulting S/N distribution for the both mass and B-mode map. 

5.2. The Impact of SkyFrame 

The most efficient procedure for reducing the B-mode systematics is "SkyFrame" . Figure |H] shows the displacement 
vector of positional difference for stars between the distortion free coordinate on the original warp and the stereographic 
projection of the celestial coordinate from the USNOA2.0 catalog at the tangent point. Although the displacement 
is small, a few pixels, there is a coherent pattern. Smoothing during construction of the map can introduce coherent 
patterns that dominate the noise from random ellipticities. These patterns can introduce the small scale B-mode 
signal. It is also possible to generate false E-mode signal from this coherent pattern. 

There are two possible explanations for the residual pattern in the original warp: (1) an anisotropic distortion 
resulting from optical aberration, and/or (2) a disturbance in the curvature of the wave front due to atmospheric 
effects. The globally rotated component results from a small rotation of the image rotator. Higher order patterns vary 
fr om the field-to - field. 

Yos hino et al.l (|2Q07f ) evaluated the anisotropy in the distortion from the corrector lens in Suprime-Cam with a 
positional accuracy of 0.4". They found that there is an ~0.8" difference in the isotropic distortion pattern between 
the right and left edges of the field-of-view. These distortions are comparable with the residual of the displacement 
from the distortion free coordinate on the original image to the p rojec ted celestia l coor dinate. 

These patterns are different for different observed fields (Figure 9(a) and Figure 9(b) ). The difference between fields 
is the same as the difference of exposure-to-exposure because all frames are simply registered to the first frame in the 
series of exposures. These issues can be the source of the patterns we find. 

A non-static factor like an atmospheric disturbance cannot explain the patterns. Ide Vries et al.1 (|2007| ) show that 
longer exposures yield progressively lower ellipticities of the PSF (as l/^/t exp ) suggesting that the coherent differences 
shown in Figure 9(a) should not be caused from atmospheric disturbance. 



6. DEFINING THE WEAK LENSING PEAKS 

At least two issues are important in defining the weak lensing peaks in the k-S/N map. The smoothing length is 
important because it reduces the shape noise. Furthermore, the smoothing in part determines the sensitivity to cluster 
detection as a function of mass and redshift. The threshold for peak finding is also important. The fraction of false 
positives is sensitive to the threshold. Here we discuss these two issues. 

6.1. Smoothing Length 

In the previous sections, we constructed weak lensing mass maps with a fixed smoothing length of 1.5 '. A smaller 
smoothing length introduces more noise in the map resulting from shape noise. A larger smoothing length (larger than 
a few arcminutes) smooths out the cluster signal; the angular diameter of typical clusters for redshifts 0.2 < z < 0.7, 
the range of sensitivity of the weak lensing maps, is se veral arcminutes 



Previous studies ad opt a range of smoothing lengths. dWittman et al.l200lUMivazaki et alJ2002al:lHetterscheidt et al.l 
20051: iWittman et all 12003 iSchirmer et alJ 120071 iGavazzi fc Soucaill l200fl iMivazaki et al.l 12001 : iMaturi et all 120071 : 



Berge et al.ll2008t iDietrich et al.ll2007l) . To examine the effect of the smoothing length, we construct E/B-mode maps 
and random maps with three different smoothing lengths: 6q — 0.5', 1.0', and 1.5 '. 

Figure [TOl shows the relation between x 2 and 9q- The x 2 of the E-mode increases with the smoothing length. The 
X 2 of the B-mode also increases. 
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Fig. 7. — Mass maps for the GTO field overlaid on Subaru images. Contours are drawn from 2<r threshold with an interval of 0.5a. 
Boxes show identified peak positions along with peak number and its significance. Masked region for brighter stars are also overlaid. 
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Fig. 8. — S/N distributions for E-mode {triangles) and B-mode {circles) maps for the GTO fields. Lower panel shows the deviation from 
the normal distribution for the B-mode map. 
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Fig. 9. — Displacement vector map of stellar positional differences between the undistorted detector coordinate (Section 12. 3H and the 
stereographic projection of the celestial coordinate from the USNOA-2.0 catalog at the tangent point for the GTO subfield. The box shows 
approximately 32' X 27' which is the field of view of Suprime-Cam. The length of the arrows indicate the positional difference in the 
detector coordinate multiplied by 1000. The displacement is complex and includes rotational components at the level of the pixel scale. 



15 

The averaged x 2 value among 100 random realizations is x 2 > 1 and it increases with the square root of the 
smoothing length 9q. This result occurs because the smoothing reduces the effective number of independent pixels. 
For small 9q the behavior of x 2 comes from the pixelization caused by dividing the smoothing kernel into coarse pixels. 
The shortest smoothing length, 9q — 0.5', corresponds to only w 1.5 pixels. The coarse sampling cannot reproduce 
the gaussian smoothing kernel well; thus the amplitude of the scaling relation may be changed and increased. This 
kind of pixelization effect is hard to describe analytically. Because the range, 0'.5 < 9q < l'-5, might be affected by 
this pixelization effect, we do not correct for the number of effective pixels. We use the scale dependent x 2 (#g) as a 
measure of the x 2 ■ 

We estimate the o ptimal smoothing length for detecting clusters as a function of mass and redshift. Following 
lHamana et al.l (|2004l) . we analytically evaluate the expected S/N for weak lensing halos with universal NFW profiles. 
Because more distant clusters have smaller angular diameter, the S/N is redshift dependent. Furthermore, the typical 
scale radius of an NFW cluster depends on its mass; thus the S/N is also mass dependent. We compute the S/N 
as a function of smoothing scale 9q. The S/N for less massive clusters (-M14 = 10 14 h~ 1 MQ ) ; M14 = 3.0 for z = 0.5, 
A/14 = 1.0 for z = 0.2) reaches a maximum at 9q f=a 1.0'; more massive clusters are detected at greater S/N with a 
smoothing length larger than 9q sa 1.0'. However, difference between 1.0' and 1.5 ' is small. At a fixed cluster mass, 
smaller smoothing lengths yield higher S/N at larger redshift. Because these effects are small, we choose 9q = 1.5' . 

6.2. Peak Definitions and Thresholds 

We use SExtractor to locate peaks in the E-mode map; we require at least 3 connected pixels (> 2a) to define a 
peak. Selecting the appropriate threshold for a clean peak catalog is s subtle issue. Figure [TT] shows the cumulative 
histogram of the highest peaks found in each noise map as a function of S/N. This figure shows the likelihood that 
significant peaks can be found accidentally. Statistics describing the highest peak for the 100 realizations are shown 
in the figure and summarized in Table |H 

We compare the cumulative peak counts with the expected standard normal distribution (green line). Figure [Til 
shows that the highest peak statistics of the noise maps are not well described by the theoretical curve. In other 
words, the highest peak statistics do not obey the ideal random Gaussian field. Thus we must derive these statistical 
quantities from Monte-Carlo simulations. 

This discrepancy between the noise maps and the ideal case is more significant for smaller smoothing lengths. Based 
on these tests, we select different thresholds for each smoothing length. We consider two kinds of threshold: the 3a 
threshold and the 99% threshold. The 3a threshold is derived by assuming the standard normal distribution for the 
highest peaks; the 99% threshold is derived from "the most extreme peaks" among the highest peaks in the entire set 
of 100 noise map. Because it is based on the real data, the 99% threshold is much more reliable than the 3a threshold 
and we use it throughout. 

We summarize the 3a and 99% thresholds in Table |4] As expected these thresholds differ for smaller smoothing 
lengths, but they are similar for 9q = 1.5. 

The numbers of detected peaks for each threshold are listed as N 3(J , N 99 % in Table |H 

TABLE 4 
S/N and Peak Statistics as a Function of the Smoothing Length 0q 





S/N values 






Max Peak 




Detected Peak 


e G 


Xr Xr c 


M 


a 


(S/N) Sa 


(S/N) 99% 


AW AW„ 


0.5 


1.94 29.9 1.01 


4.07 


0.25 


4.82 


5.05 


2 


1.0 


15.1 175 1.05 


3.83 


0.22 


4.48 


4.64 


4 4 


1.5 


54.4 545 1.11 


3.72 


0.30 


4.62 


4.57 


4 4 



a "S/N values" shows statistical properties for the S/N histogram, the merit functions for E/B-mode, and the standard deviation of pixel 
values. "Max Peak" shows statistical properties for the highest peak value (S/N)3 a and the threshold given the most extreme peak among 
the highest peaks in the entire set of 100 noise map {{S/N) 9 g%). "Detected Peak" show the number of detected peaks corresponding to 
the adopted threshold. 

Next we consider the peaks in the E-mode map. With the smallest smoothing, there are no peaks. The shape 
noise of intrinsic ellipticities raises the threshold. The situation is different for Bq — 1.0 and 1.5 because the larger 
smoothing involves more sources. The lensing signal becomes more significant for larger smoothing lengths. 

As a quantitative representation of the significance of the weak lensing signal, we evaluate X e /Xb- For the B-mode, 
X 2 should be small because there should be no signal; for the E-mode, x 2 should be larger reflecting the lensing signal. 

The actual values of Xe/Xb f° r @G — 1-0 and 1.5 are 11.6 and 10.0 for the GTO field. These similar values imply 
there is no strong preference for smoothing lengths between 1.0 and 1.5 in our data. 

7. THE DLS KAPPA MAP 

The procedures of Sections 15.11 generally apply to both the GTO and DLS field maps. Here we discuss the few issues 
that are particular to the DLS map. 
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Fig. 11. — Cumulative distribution of the highest peak S/N values for 100 realizations of noise maps for the GTO field. The solid line 
shows the expectation from the normal distribution. 

For the DLS map, the x 2 f° r the original procedure is 134; for the new one it is 49.1 ( Figure IT21 . The new value 
although reduced from the original one is not within the range for the random realizations (32 is the maximum value). 
The S/N distribution for the B-mode has a slight excess on the negative side. 

We use i' band imaging in the DLS field rather than the Re band. We cannot determine the PSF anisotropy near 
the edge of the field because the star/galaxy separation is poor. Thus we trim objects located at the outer edge (> 15' 
of the center of field). In contrast we do not trim the edges of the GTO field. 

Even after this trimming, there is still a small excess in the S/N distribution for the B-mode on the negative side. 
However, the plus side seems to be well matched to the standard normal distribution. If we mask the three prominent 
negative excesses, the B-mode S/N distribution is in good agreement with the standard normal distribution. Thus the 
global correction may be sufficient and the residual local excess may be related to the most prominent peaks in the 
E-mode map. Figure Q2] shows the k-S/N map for the DLS field. Here we present the resulting k-S/N maps for the 
DLS field (Figure [£3]) 

8. DISTRIBUTION OF PEAKS IN NOISE AND B-MODE MAPS 

It is useful to have an analytical expression for the distribution of the highest peaks in order to test both the noise 
and B-mode maps. These expressions can then be used to estimate the fraction of false peaks in the k-S/N maps. 

We derive a probability density function that enables us to evaluate the probability of finding the highest peak in the 
noise map as a spurious peak in the k-S/N map. We evaluate this probability as a function of the S/N. The probability 
density function that a random value x obeying the normal distribution becomes the maximum value among n samples 
is: 



P n (x)dx 



exp 



4 G- 



V2j 



n-\ 



dx 



(8) 



where erfc(x) is the complementary error function: erfc(x) = Ij^pH \ e * dt. The first part of the product gives the 
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Fig. 12.— Same as Figure|8]but for the DLS field. 
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probability of obtaining a random value of x, and the latter part describes the probability that the remaining n — 1 
random values are < x. 

However, our map creation procedure is more complex because it applies gaussian smoothing which correlates pixel 
values spatially. We have to check the effect of the spatial correlation prior to evaluating the probability. We apply the 
smoothing function on the k map in the form exp (— (8/6q) 2 ) where 8q/\/2 is the standard deviation of the smoothing 
kernel. Using the standard deviation as a measure of the correlation length, we reduce the total number of pixels in the 
mass map by a factor of (9q/^/2) 2 . For example, in the GTO field where the total number of pixels is 300 x 234 with 
a 0.33' / pixel scale, the effective number of pixels as a result of the smoothing becomes 300 x 234 x (0.33 x -\/2/#g) 2 - 
After adopting 8q = 1.5', the effective number of pixels is 6,795. 

Figure Q3] shows both the analytic curve and the Monte Carlo results based on 10000 random realizations of the 
noise map. The analytic expression is an excellent representation of the Monte Carlo results. 

In summary, Figure [T4l shows that the smoothing procedure only reduces the number of independent pixels. Using 
this effective number in the probability function, we can evaluate the confidence level corresponding to a given S/N 
threshold. For example, a threshold of S/N = 4.66 corresponds to a 99.7% confidence level for the 2.12 deg 2 B-mode 
free massmap based on 1.5' smoothing. Similarly, S/N = 4.50 and 5.80 are the same 99.7% level confidence levels for 
a 1 and 1000 deg 2 degree survey, respectively. 

Using this same formalism, we can derive an expression for the fraction of spurious peaks in the k-S/N map above 
a given threshold. The probability of obtaining the (to + l)-th value in the descending set of ordered random values 
as a function of S/N is: 

,« ( ^.fc^^« P (4)(^ fc (^))"-"'- , ( 1 4„ fc (^))"' dI „, 



where n C m = n\/{n — m)\m\. With this general probability function, we can evaluate the expectation value for the 
number of peaks found in the noise field with a given S/N threshold: 

/•oo n — 1 

N = £4"° (*)<** ( 10 ) 

J(S/N)th m=0 

According to this formula, we should find on average 2.7 and 25.3 spurious peaks exceeding (S/N)^ = 4.5 and 4.0, 
respectively, for a 1000 deg 2 survey with 1.5' Gaussian smoothing. 

In conclusion, the analytical expression we derive is useful for future error analysis of large surveys. If we place a S/N 
= 4.5 threshold on the B-modc-frcc kappa (S/N) map, we expect these high significant peaks to be uncontaminated 
at the 99.7% confidence level for a 1 deg 2 field; only 2.7 spurious peaks should be found on average for a 1000 deg 2 
survey. Figure [15] shows the S/N threshold as a function of survey area. 

9. THE NATURE OF THE WEAK LENSING PEAKS 

To test the reliability of the peak catalogs for the GTO (Figure [7]) and DLS fields (Figure fT3"|) . we next search for 
corresponding systems of galaxies. Weak lensing only traces the surface mass density. Thus our peaks may be single 
isolated clusters or superpositions of less massive groups alon g the line of sight. Here we test the peaks against dense 
foreground redshift surveys following Gel ler et"aLl ()2005D and lKurtz et al.l (|2012t ). 



9.1. Redshift Surveys 

Dense redshift samples are available for both the GTO field (|Kurtz et al.ll2012l ) and t he DLS field dGeller et al.ll2005l 
120101 : [Hwang et aT]|2012j ). The redshift surveys were carried out with the Hectospec (jFabricant et al.lll998l 120051 ). a 
300-fiber robotic instrument mounted on the MMT. The spectra cover the wavelength range 3, 500-10, 000A with a 
resolution of ~ 6 A. The typical error in an individual redshift (normalized by (1 + z)) is 27 km/s for emission line 
objects and 37 km/s for absorption line objects. The surveys include redshifts for galaxies with R < 20.6 in the 4 
deg 2 DLS field and redshifts for galaxies with r pot ro < 21.3, r - i > 0.4 and g — r > 1.0 in the GTO 2deg 2 field. With 
these samples, we can identify massive clusters of galaxies with redshift z < 0.6. The surface number densities of the 
spectroscopic samples are 0.6 and 1.1 arcmin -2 for the GTO and the DLS fields, respectively. 

To identify potential galaxy systems associated with the detected weak lensing peaks, we examine the redshift 
distributions in cones of 3' and 6' radius centered on the location of each of the weak lensing peaks. 

We examine peaks in the redshift histograms. We define the significance of each peak as the difference between the 
number of galaxies N and mean number of galaxies (N) expected in the cone at that redshift divided by the standard 
deviation a: (N — (N))/a. For neighbors satisfying |Az| < 0.004(1 + z), we calculate a velocity dispersion from the 
square root of the unbiased variance. We use the bootstrap to measure the error in the velocity dispersion. 

Tables [5] and [5] list the weak lensing peaks for the GTO and DLS fields, respectively. The Tables list the lensing 
peak rank and significance, the celestial coordinates of the peaks, the rest frame velocity dispersion and its error, the 
number of peak members and its significance, and the mean redshift. The Table also include additional information 
from the literature. 
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Fig. 14. — S/N value distribution for maximum peaks among 10,000 noise map realizations for a field covering 2.12 deg 2 (histogram) and 
an analytical model (curve). The effective pixel number is the number of actual pixels divided by the smoothing length in pixel units. 



9.2. The GTO Field 



We list peaks in Table [5] that have S/N > 3.7 as in [Mivazaki et all ()2007l ). However, we concentrate on peaks with 
S/N> 4.56 because our analysis of the noise maps shows that this higher threshold yields a more robust catalog. 

The most si gnificant peak, 0, has an obvious counterpart in the redshift histogram at a mean redshift z = 0.419 
(Figure 16(a) I. Figure 16(b) shows the spatial distribution of cluster members within Az = 0.004 of the mean and 
within &' ot the weak lensmg peak. It is clear that these galaxies are concentrated around th e lensing peak. T here 
is also extended x-ray emission associated with this system; the SDSS red sequence survey (|Wen et al.l I2009D also 
identifies it. 

The lensing peak is slightly offset from the galaxy distribution. Assuming the brightest cluster member is the BCG, 
its offset is ~ 1', comparable with smoothing length. 

Figure [T71 corresponding to peak 1, shows a concentration of bright galaxies, but the redshift histogram does not 
show an obvious pea k. This discrepancy may result fr o m un der-sampling at this location in the redshift survey (see 
peak 7 in Figure 2 of lKurtz et al.l (|2012[ )). iGunn et al.l ()1986l ) reported that this concentration of galaxies is a cluster 
with z = 0.2830. The redshift histogram does show some galaxies around this redshift. 

Peak 2 has a significance of 4.92 (Figure IT8"|). This peak has a clear c ounter part in the redshift histogram at a mean 
z = 0.540. X-ray observations also show extended emission. Figure 18(b) shows the distribution of the suspected 
cluster members on the sky. The lensing peak is shifted by several aremmutes from the peak of galaxy distribution. 
The reason for the shift is not clear, but we suspect that because the peak is near the edge of the imaging field, there 
may be distortion in the k-S/N map. This peak may also be affected by a neighboring bright star. 

Peak 3 has a significance of 4.76 and, remarkably, corresponds to a peak in the redshift survey at z = 0.73 (Figure 
TS|) . Although the spectroscopic sample is small, there is an excess in the galaxy count within 3' of the peak. 

Above the S/N = 4.56 threshold, all of the peaks (with the possible exception of peak 1) correspond to systems of 
galaxies. We suspect that the failure to detect peak 1 in the redshift survey results from undersampling. The offsets 
in position between the lensing peaks and the lenses are consistent with the smoothing length except in a case where 
edge effects may be important. 

Below the S/N = 4.56 threshold many fewer peaks correspond to systems of galaxies as we expect based on analysis 
of the noise maps. Table [5] lists possible systems associated with S/N> 3.7 peaks for comparison with previous work. 



9.3. The DLS Field 
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Fig. 15. — S/N threshold dependence on survey area. 
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(a) Redshift distribution for Peak GTO 00. The filled histogram (b) Close-up view for Peak GTO 00. The image is 6' X 6'. Con- 
shows objects within a cone of 3' radius centered on the weak tours show the weak lensing peak. Small yellow circles identify 
lcnsing peak; the open histogram shows objects within a 6' cone, galaxies in the redshift peak centered at a mean z = 0.416. 
Bins in redshift are 0.002(1 + z) wide. 



Fig. 16.— Peak GTO00 



TABLE 5 
Peak list for the GTO field. 



Rank 


V 


RA2000 


DEC2000 


o- Pi6 [km/s] 


N 6 


Z6 


o- p ,3[km/s] 


N 3 


Z3 


NED 


X-ray 


distance 





5.95 


16:02:50.826 


+43:35:51.34 


992.0 + 94.5 


40(8. 60-) 


0.415(#1) 


1014.1 + 113.2 


24(11. 8<r) 


0.416(#1) 


0.420^, XrayS^ 


9.10 + 1.32* 


1.156' 


1 


5.07 


16:07:58.418 


+43:38:23.99 


.a 






.a 






0.2830 c 






2 


4.94 


16:03:01.380 


+42:46:32.00 


655.8 ± 94.8 


29(12.0<t) 


0.540(#0) 


728.3 + 127.7 


18(13. 9cr) 


0.540(#0) 


0. 556^,0. 539 V_,XrayS^ 


20.0±1.5 f 


1.931' 


3 


4.76 


16:05:24.866 


+42:45:34.17 


497.9 ± 136.5 


10(16.1<t) 


0.731 


562.8 + 151.1 


8(26.2cr) 


0.731 








4 


4.38 


16:05:39.982 


+43:22:13.70 




















5 


4.29 


16:04:37.615 


+43:27:34.65 














XrayS^, 






6 


4.19 


16:04:41.282 


+42:37:54.73 




















7 


4.17 


16:07:44.303 


+42:37:45.58 




















8 


4.12 


16:04:10.464 


+42:39:14.47 


411.9 + 74.9 


9(18. lcr) 


0.673(#3) 


508.6 ± 102.4 


6(28.1cr) 


0.673(#3) 








9 


3.85 


16:05:15.727 


+42:39:14.39 




















10 


3.70 


16:01:46.456 


+42:56:06.40 


421.9 + 59.1 


15(4.5cr) 


0.256 


366.1 + 59.7 


10(7.2cr) 


0.287 









The column "Rank" is sorted by the lensing significance v. a v shows a measured velocity dispersion, a number of galaxies on the cone and its mean redshift for 6 and 3 arcmin 
cones, respectively. NED: objects classified as "Galaxy Clusters" or "X-ray sources (XrayS)" in NED within 3 arcmin, The column "X-ray" and distance shows a flux of X-ray in the 
literature and a positional difference between detected peak and X-ray source. 

a This peak is not well sampled at this location in the redshift survey 

b WHO: HWen et a l. 2009, 2010), 

c GHO: (Gurm ct al. 1986), 

d lWGA: ftp: //cdsarc.u-strasbg.fr /cats/IX/30/ReadMe 

°KLG2009: IjKocevsk i ct al. 2009) 

f X-ray flux is converted from count rate in 0.24-2.0 keV with a constant correction factor of 1.5 X 10 — 11 erg/cm 2 /sec in the unit of 10 _14 erg/cm 2 /scc (White ct al. 2000) . 
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(a) Same as Figure |16(a)] but for Peak GTO01. 
Fig. 17.— Peak GTO01 
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(b) Close up view for Peak GTO01. 
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(a) Same as Figure |16(a)] but for Peak GTO02. 



(b) Sa me as Figure |16(b)| The mean redshift of the peak 
Figure |18(a)| corresponding to the weak lensing system is z 
0.540 



Fig. 18.— Peak GTO02 



In th e DLS field (Figure [j~3j). the most significant peak (0) has S/N = 6.18 (Figure l20|) . The redshift histogram 
(Figure 20(b) I shows two different concentrations centered at z — 0.320 and z — 0.537. This peak is apparently a 
superposition of clusters along the line of sight. So far, this peak is not covered by X-ray observations. 

Peak 1 (Figure |2"TJ) has S/N = 4.98 but it unfortunately coincides with bright stars. The center of the weak lensing 
peak is masked by these bright stars. The redshift histogram shows a peak at z — 0.531, but the velocity dispersion 
for the 3' cone, a ~ 386 km s _1 is too low to account for the lensing peak: the 636 km s" 1 dispersion for the 6' 
cone is sufficient to account for the detection. The coincident star compromises accurate measurement of the velocity 
dispersion. 

Two other previous studies (lensing and X-ray) (|Abate et al.ll2009t iKubo et"aTll2009D reported that there is a cluster 
at z = 0.53. X-ray observations show extended emission suggesting that the velocity dispersion may be underestimated. 
The redshift distribution also shows a peak at at z = 0.183 Thus the system at z = 0.53 may be boosted by the 
foreground large-scale structure. 
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(a) Same as Figure |16(a)] but for Peak GTO03. 
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(b) Same as Figur e [J6(b)| The significant peak in the redshift 
histogram (Figure |ljj(aj[ | corresponding to the lensing peak is 
2 = 0.731. 



Fig. 19.— Peak GTO03 
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(a) Same as Figurc |16(a")] but for Peak DLS00 
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(b) Same as Figure [J6(b)| This weak lensing peak corresponds 
to a superposition ot structures in the redshfit survey at red- 
shifts z = 0.319 (yellow circles) and z = 0.535 (skyblue circles). 



Fig. 20.— Peak DLS00 

The third significant Pea k 02 (F igure [22]) has S/N = 4.56 and coincides with five galaxies around z = 0.65 in the 
spectroscopic data (Figure 22(b) I even though such high z galaxies are rare in the redshift survey. There is also 
extended x-ray emission associated with this peak. 

In the DLS field, the results are less clear than in the GTO field. One of the peaks (2) is a cluster and one (0) is a 
clear superposition of clusters. The interpretation of peak 1 is unclear, but we believe the peak is not noise. In Table 
[6] we also list lower significance peaks and possibly associated systems of galaxies for comparison with previous work. 

10. COMPARISON WITH EARLIER RESULTS 

The new maps of the GTO field and our 1 deg 2 portion of the DLS field contain seven weak lensing peaks above 
our threshold; these peaks are probably all associated either with a single system of galaxies or with a superposition 
along the line-of-sight. In other words, the analysis suggests that the lensing signal above our threshold is a reasonably 



TABLE 6 
Peak list for the DLS field. 



C: 



Rank 


V 


RA2000 


DEC2000 


o- p , 6 [km/s] 


N 6 


Z6 


°>.3[km/s] 


N 3 


zz 


NED 


X-ray 


distance 





6.18 


09:16:52.817 


+30:00:28.44 


497 ± 50 
486 ± 75 


47(7.8cr) 
25(7.0o) 


0.318 
0.535 


468 ± 61 
446 ± 79 


23(9. 3cr) 
18(13.7cr) 


0.319 
0.535 








1 


4.98 


09:15:54.449 


+29:37:08.25 


636 ± 133 a 


26(7.2cr) 


0.531 


386 ± 64 a 


12(9.2o-) 


0.531 


0.53 bc 


4.99 + 1.32(2.97)1 


0.797' 


2 


4.65 


09:16:04.967 


+30:28:08.30 


524 ± 145 


5(30.2o) 


0.647 


524 ± 147 


5(54. lo) 


0.647 




9.52 + 2.24(4.60) 
1.75 + 1.21(0.65)1 


0.786' 
1.677' 


3 


4.46 


09:15:54.292 


+30:03:08.24 




















4 


4.24 


09:17:12.960 


+30:18:28.09 


811 ± 124 


19(5. 3cr) 


0.535 


859 ± 162 


12(9. lo) 


0.535 




29.1 ± 10.4(5.88)1 


0.961' 


5 


3.86 


09:16:20.559 


+29:15:48.60 


1026 ± 188 
818 ± 99 


13(11.2<r) 

25(7.0cr) 


0.365 
0.536 


486 ± 187 
968 ± 169 


8(17.7o) 
12(9. lo) 


0.367 
0.536 









Same as Table[5]but Peak list for the DLS field. NED: objects classified as "Galaxy Clusters" or "X-ray sources (XrayS)" in NED within 3 arcmin 
a This peak coincides a foreground bright star and then prevent us to measure accurate velocity dispersion. 
b AWM2009 (Abate ct al. 2009) 

C KKD2009 UKubo et al.H2009f) . 

d X-ray flux in 0.2keV-12keV(0.2keV-2keV) band for extended soruce (ext> arcsec) in the unit of 10 — 14 erg/cm 2 /s (Xmm-Ncwton Survey Science Centre 2010) 
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(a) Same as Figure [16(a) | but for Peak DLS01 



Fig. 21.— Peak DLS01 




(a) Same as Figure [16(a)] but for Peak DLS02 
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(b) Sa me as Figure |16(b)| The mean redshfit of the peak ir 
Figure |21(a)| corresponding to the weak lensing system is z = 
0.531. 
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(b) Sa me as Figure |16(b)| The mean redshfit of the peak 
Figure |22(a)| corresponding to the weak lensing system is z 
0.647. 



Fig. 22.— Peak DLS02 



robust measure of features of the large-scale structure of the universe. 

Next we compare the results for the lensing peaks with previous work. Th e issues relevant fo r the GTO and DLS 
field are slightly different. In the GTO case, we analyze the same data as Kurtz et al.l (|2012f) . but we apply new 
procedures to reduce systematic error. In the DLS subfield, we apply essentially the same reduction procedure as the 
GTO field. Ou r Suba ru DLS data are deeper and are taken in better conditions than the original data analyzed by 
IWittman et"all (|2006ft 



For the GTO field. IKurtz et al.l (J2012D find two peaks with S/N > 4.5; we also find these peaks at S/N > 4.56 along 
with two additional ones. This difference occurs because the suppression of systematic errors changes the local peak 
values in the re-S/N map. There is no change in the average sensitivity of the map. 

The difference in the peak height distribution for the B-mode map of IKurtz et al.l ([2012D and our map demonstrates 
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again that the systematic error suppression is effective. iKurtz et al.l (|2012l) found 4 peaks with S/N> 3.7 in their 
B-mode map. We find only 2 peaks with this amplitude in our new B-mode map (S/N=3.90, 3.87). The highest 
B-mode peak now has S/N = 3.90 rather than 4.23. 

Our data cover about 25% of the original DLS field F2. iWittman et al.1 (|2006[ ) identified one peak and lKubo et al.1 
(2009) detected 3 peaks in this region of F2. All of these peaks have S/N< 3.9, substantially below our revised 
threshold. 

The original DLS F2 data were taken on the Kitt Peak 4-meter in the R-band image with PSF < 0.9"; final co-added 
images reached a depth of R ~ 26 and the source galaxy density is, ~ 20 arcmin~ 2 . The Subaru source density is 
~ 30 arcrnin. - to roughly the same magnitude limit. The lower source galaxy density may result from the poorer 
seeing at Kitt Peak. In this case, the Subaru map is more sensitive mainly as a result of the increase source density. 
We find three peaks above the threshold S/N = 4.56; one of these is a cluster at z = 0.64. 

Obviously, at S/N below our revised threshold, there are peaks that correspond to systems of galaxies (Tables [5j EJ . 
As we would expect, the fraction of "real" peaks declines with decreasing S/N. 

11. CONCLUSION 

Weak lensing is a powerful tool of modern cosmology. Its application to the detection of clusters of galaxies is limited 
to some extent by subtle systematic errors in the construction of the k-S/N maps. Careful analysis of the B-mode maps 
where there should be no signal provides a route to reduction of these systematic errors. We use a set of noise maps 
along with the B-mode maps for two fields observed with the Subaru telescope to test procedures for the reduction of 
systematics. 

Our tests of the B-mode maps against noise maps demonstrate that the following image processing procedures 
significantly reduce systematic error in the final re-S/N map: 

• Registration of positions onto the absolute celestial coordinate grid, 

• Simultaneous application of all image transformations to avoid loss of information in interpolation, and 

• Application of a smoothing cutoff at 10' to remove excess large-scale power. 

By comparing the S/N distribution for noise and B-mode maps, we demonstrate that these revised procedures signif- 
icantly reduce the systematics that otherwise appear in the B-mode map. We then take advantage of these revised 
procedures to show that weak lensing cluster detection is insensitive to the smoothing length for the k-S/N map in 
the range 1.0' to 1.5'. Also on the basis of this analysis, we suggest use of the 99% threshold derived from the "most 
extreme peaks" among the highest peaks of the noise maps rather than the more standard 3<r threshold. We provide 
an analytic expression for the distribution of these highest peaks that can be used to estimate the fraction of false 
peaks in the re-S/N maps as a function of the detection threshold. 

We test our systematic error reduction procedure and our high peak analysis by analyzing the GTO field and a portion 
of the DLS F2 field; together these fields cover 2.52deg 2 . A threshold of 4.56 corresponds to the 99.5% confidence 
level from the highest peak analysis for our 1.5' smoothing length. We expect these peaks to be uncontaminated by 
false peaks over this total area. The data show seven peaks above this threshold; all of these peaks correspond to 
galaxy overdensities. Six of these peaks correspond to clusters of galaxies; one is a superposition of systems. Five 
of these peaks are well-sampled by dense foreground redshift surveys providing reliable estimates of the rest frame 
line-of-sight velocity dispersions of the clusters. One peak in the GTO field is in region sparsely sampled by the redshift 
survey, but a concentration of galaxies is evident in the imaging data. For one peak in the DLS field, foreground stars 
probably lead to an underestimate of the velocity dispersion from the spectroscopic data, but the peak is associated 
with extended x-ray emission. Taken together, these results substantiate the efficacy of our revised procedures. 

The systematic error reduction procedures we apply are general and can be applied to future large-area weak lensing 
surveys. Our high peak analysis suggests that with a S/N threshold of 4.5, there should be only 2.7 spurious weak 
lensing peaks even in an area of 1000 deg 2 where we expect ~ 2000 peaks based on our Subaru fields. 
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APPENDIX 

A. MASKS 

Here we present details of masks describ ed in Secti on [3] Tab le KT1 shows the mask for brights stars based on the 
external catalog USNO-A2.0 ()Monetlll998t ). Table 1X21 and Table POl show the coordinates for the additional manually 
constructed masks. 
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TABLE Al 
Masking Radius Around Bright Stars 



Magnitude Range 


Mask Radius (arcscc) 


R < 9.0 


210 


9.0 < R< 10.0 


140 


10.0 < R < 12.0 


70 


12.0 < R < 14.0 


35 


14.0 < R < 16.0 


18 


16.0 < R < 18.0 


6 



TABLE A2 
Manually Masked Regions around Bright Stars in the GTO Field. 



RA2000 


DEC2000 


radius (arcscc) 


Origin 


240.11876 


42.9889 


270 


Filter reflection 


240.327 


42.7408 


270 


Filter reflection 


240.44033 


42.842131 


220 


Filter reflection 


240.45127 


43.480879 


270 


Filter reflection 


240.48183 


43.629343 


41.4 


Incorrect position 


240.77826 


42.616311 


270 


Filter reflection 


241.06437 


43.414799 


99 


Dewar window reflection 


241.16981 


43.42311 


36 


Incorrect position 


241.26726 


43.561169 


99 


Dewar window reflection 


241.3551 


43.711462 


270 


Filter reflection 



TABLE A3 

Manually Masked Regions Around Bright Stars in the DLS Field. 



RA2000 


DEC2000 


radius (arcscc) 


Origin 


138.888195 


29.792573 


120 


Dewar window reflection 


139.343445 


29.452575 


120 


Dewar window reflection 


139.100387 


30.040381 


120 


Dewar window reflection 


139.254212 


30.374562 


120 


Dewar window reflection 


139.412442 


30.320542 


120 


Dewar window reflection 



